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Abstract 

Background: Sugars modulate expression of hundreds of genes in plants. Previous studies on sugar signaling, using intact 
plants or plant tissues, were hampered by tissue heterogeneity, uneven sugar transport and/or inter-conversions of the 
applied sugars. This, in turn, could obscure the identity of a specific sugar that acts as a signal affecting expression of given 
gene in a given tissue or cell-type. 

Methodology/Principal Findings: To bypass those biases, we have developed a novel biological system, based on stem- 
cell-like Arabidopsis suspension culture. The cells were grown in a hormone-free medium and were sustained on xylose as 
the only carbon source. Using functional genomics we have identified 290 sugar responsive genes, responding rapidly 
(within 1 h) and specifically to low concentration (1 mM) of glucose, fructose and/or sucrose. For selected genes, the true 
nature of the signaling sugar molecules and sites of sugar perception were further clarified using non-metabolizable sugar 
analogues. Using both transgenic and wild-type A. thaliana seedlings, it was shown that the expression of selected sugar- 
responsive genes was not restricted to a specific tissue or cell type and responded to photoperiod-related changes in sugar 
availability. This suggested that sugar-responsiveness of genes identified in the cell culture system was not biased toward 
heterotrophic background and resembled that in whole plants. 

Conclusions: Altogether, our research strategy, using a combination of cell culture and whole plants, has provided an 
unequivocal evidence for the identity of sugar-responsive genes and the identity of the sugar signaling molecules, 
independently from their inter-conversions or use for energy metabolism. 
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Introduction 

Daily changes in the environment of green plants result in 
constant fluctuation in the contents of endogenous sugars, which - 
besides being the initial carbon backbone for biological molecules - 
represent the main energy source in both the autotrophic and 
heterotrophic tissues. Primary sugars in plants include glucose 
(Glc), fructose (Fru) and sucrose (Sue), which derive from 
photosynthesis and storage reserves. Their fine-tuned distribution 
from the carbon-fixing and reserve cells to the rest of the organism 
is essential for the optimal plant growth and development [1]. 
Besides their role as energy resource, these sugars also act as 
signaling molecules, modulating the expression of the so called 
"sugar-responsive" genes. In plants, there are hundreds of such 
genes, involved both in the cellular adjustment to nutrient 
availability and controling resource distribution between tissues 
and organs . Studies on sugar-responsive genes used a variety 
of biological experimental systems, ranging from mesophyll 
protoplasts [5] and cell culture [6], to exposing the seedlings to 
a given sugar-containing media [7,8], feeding detached leaves with 



a given sugar ([9,10] and studying gene expression in whole plants 
under conditions of physiological depletion of endogenous sugars 
(e.g. end of the night) [4] or in mutants with impaired endogenous 
sugar levels [11]. Based on some of those studies, independent 
signaling pathways have been postulated for Glc, Fru and Sue 
[5,8,10]. 

When feeding specific sugars to intact plants or excised leaves, it 
is frequently not clear how fast and evenly the taken up 
carbohydrates are transported within the tissue and whether a 
given sugar is converted into related sugar species, which in turn 
might act as a distinct signal inducing changes in gene expression 
([12]. Cells lying on the exterior of the treated plant face a higher 
sugar concentration than cells on the inside of the plant, leading to 
gradients in gene expression from cell layer to cell layer. Also, a 
sample derived from a whole seedling or leaf rosette always 
consists of cells in different developmental stages, possessing a 
variety of specialized functions and distinct metabolism. This 
complicates the distinction between direct and indirect conse- 
quences of sugar treatment. A homogenous response would be 
expected with cell culture, but those cells usually require an 
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addition of both sugars and hormones to sustain their division and 
growth, which could affect sugar perception and metabolism. 
Studies using a mesophyll protoplast transient expression system 
[5] pardy overcome this difficulty; however, these photosynthetic 
cells - devoid of cell walls and composed of multiple cell types (e.g. 
spongy and palisade mesophyll) - represent a system which cannot 
be reliably extended to whole plants when studying global 
expression. 

The choice of proper experimental system to study signaling, 
whether by sugars or other compounds, represents a serious 
challenge in all fields of biology. In the medical field, stem-cell 
research has allowed for pioneering studies by using a homoge- 
nous cell population (without interference from other cell types 
and/or morphogenic compounds), which could be induced to 
differentiate into specialized cell-types, e.g. neuron [13] or vascular 
precursor ([14]. Stem-cell studies in plants have similarly enabled 
breakthroughs in understanding key compounds controlling the 
cell status from active proliferation to differentiation [15,16]. In 
the present study, in order to bypass the previously-mentioned 
biases inherent to studies on sugar signaling, we have developed a 
novel system, based on the use of a stem-cell-like suspension 
culture of Arabidopsis thaliana which could be grown without 
addition of any hormones [17]. Using functional genomics in this 
biological system, we have identified and characterized genes that 
rapidly respond to low sugar concentration and have unravelled 
true nature of the signaling sugar species and their perception sites. 
The sugar responsiveness of selected genes was then verified in 
planta using Arabidopsis seedlings. The findings, combining the use 
of habituated cell culture and whole plants, shed new light on 
sugar signaling in plants. 

Materials and Methods 

Cell culture maintenance and growth characteristics 

Habituated stem-cell-like suspension culture used in this study 
was derived from roots of A. thaliana [17]. The cells were grown on 
an orbital shaker in the dark in full-strength MS liquid media. 
Commercially available MS salts with vitamins (Duchefa, 
Haarlem, Netherlands) were prepared and autoclaved prior to 
the addition of the carbon source. Filter sterilized stock solutions of 
Sue and the pentose D-(+)-Xyl (Sigma-Aldrich) were added to the 
liquid media to reach a final concentration of 87 mM. The cells 
were subcultured every 7 d through harvesting cells by centrifu- 
gation for 2 min at 200xg, measuring the packed cell volume 
(PCV), washing the cells with fresh media and subdividing the 
culture in additional fresh media with a ratio 1:2 (Xyl-grown) or 
1 :8 (Suc-grown). To estimate the effect of the C-source on the 
growth, A. thaliana cells, conditioned for growth on either Sue- or 
Xyl-containing MS media, were subcultured into the respective 
Sue- or Xyl-containing media and incubated for 7 days in the 
dark. Every day, samples were taken and the PCV of each sample 
was estimated after centrifugation for 2 min at 200xg. Subse- 
quently, the media was removed through vacuum filtration and 
the cells were washed twice using full strength MS-media without 
C-source. The dry cell pellet was immediately frozen in liquid N 2 
and stored at — 80°C for further analysis. 

Cell culture induction conditions 

All cell culture induction assays to estimate the effect of 
exogenously supplied sugar or sugar analogues were done using 7- 
d-old Xyl-grown cells, which were not sub-cultured or transferred 
to fresh media prior to the experiment. Homogeneous cell 
solutions were equally divided onto conical flasks prior to the 
treatment. Induction of cells with 0, 1 and 5 mM filter-sterilized 



Sue was conducted for 0, 5, 15 and 30 min as well as 1 and 3 h on 
an orbital shaker in the dark. Samples for the microarray study 
were collected after a 1 h incubation of Xyl-grown cells with 
1 mM filter sterilized Xyl, Sue, Glc and Fru, respectively. The 
effect of non-metabolizable sugar analogues on gene expression 
was tested through the application of 1 mM filter-sterilized Xyl, 
Sue, D-Glc, Fru, turanose, palatinose, 2-deoxyglucose, L-Glc and 
mannitol to 7-d-old Xyl-grown cells, respectively. After treatment, 
the culture media was removed through vacuum filtration. The 
cell pellet was washed 2 times with full strength MS-media without 
sugar and, after drying, immediately snap-frozen and ground in 
liquid N 2 . The cell powder was stored at -80°C until further use. 

Plant material 

Seedlings of Arabidopsis thaliana col-0 wt, pgml mutant and 
promoter::GUS lines for the genes bZIP63, MRS2.11, URH1 (Col-0 
background) and TOR (VVassilijewska background) were grown in 
vitro on vertical plates containing a half-strength MS agar without 
sugar supplement. The harvest was carried out at the develop- 
mental stage of the appearance of the second pair of real leaves, 
which occurred after 13 days of growth in a long day growth 
chamber (16 h light, 8 h dark; 150 |Xmol m 2 s 2 ). Seedlings used 
for RNA extraction were harvested at the end of the light period 
and at the end of a 6 h extended night, by separating roots from 
shoots and immediately snap-freezing the material in liquid N 2 . 
Seedlings used for GUS-analysis were harvested at the end of the 
light period and directiy transferred into GUS containing staining 
solution. 

RNA extraction and cDNA synthesis 

Total RNA from samples after Sue-treatments and sample 
preparation for the microarray was isolated by applying TRI 
Reagent (Sigma-Aldrich) and chloroform to frozen and ground 
plant material. Phase separation was reached by centrifugation at 
full speed for 15 min. The aqueous fraction was mixed with half a 
volume isopropanol and the nucleic acids were pelleted by a 
1 0 min centrifugation at full speed. The pellet was washed 2 times 
with 70% ethanol, air-dried and subsequently taken up in sterile 
water. The total RNA samples were DNase-treated using RQ1 
RNase-Free DNase (Promega Biotech AB) and RNA was re- 
extracted as described above. Total RNA from samples after sugar 
analogue treatments as well as seedling material was isolated by 
using the E.Z.N.A Plant RNA kit from OMEGA BioTek (R6827- 
01). After elution of the RNA from the column, possible DNA 
contaminations were removed through a DNase treatment using 
the Ambion DNA^/ree DNase treatment and removal kit 
(AM1906). cDNA was generated using Superscript II Reverse 
Transcriptase (Invitrogen) and Ohgo(dT) 18 primers (Fermentas 
Life Science) following the manufacturers protocol. 1 fig total 
RNA was used in each reaction. 

Gene expression analysis 

Quantitative Real-time PCR was conducted using Roche 
LightCycler 480 Real Time PCR system and Roche LightCycler 
480 SYBR Green I Master following the manufacturer's protocol. 
Sequences of RT-PCR primers for all tested genes are listed in 
Table SI. TUB5 was used as endogenous reference gene in 
samples derived from A. thaliana cell culture after analysis of the 
microarray data showed no changes in the expression in response 
to sugar treatment. In samples derived from A. thaliana seedlings, 
PP2A was used as endogenous reference gene. The shown data 
represent means ±SD of the normalized expression in relation to 
the expression of the respective reference gene. To assess the 
statistical significance of differences in gene expression between the 
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various treatments, a Student's Hest (two-tailed distribution of two 
samples with unequal variance) was applied. 

The analysis of the transcriptome, using the Affymetrix 
GeneChip ATH1 genome array on 3 biological replicates per 
treatment was done by the SciLife Lab (Uppsala, Sweden). The 
resulting signal values were transformed into expression data and 
normalized using the RMA normalization algorithm. Subsequent- 
ly, the ratios of the gene expression of each gene in Sue-, Glc-, Fru- 
to Xyl-treated cells were calculated, respectively. Finally p-values 
were calculated using Student's /-test. Genes showing p-values < 
0.002 were selected as differentially expressed and sorted into non- 
redundant lists with genes responsive to Suc/Glc/Fru, Suc/Fru, 
Suc/Glc, Fru/Glc, Sue, Glc or Fru, respectively. Hierarchical 
clustering based on the Pearson correlation of the expression 
profile of the selected genes was performed using the Hierarchical 
Clustering Explorer v.3.5 software. The list with selected 
differentially expressed genes was cross-referenced with lists of 
sugar-responsive genes deriving from other global studies listed in 
Excel file (Dataset SI). In addition the expression of the selected 
differentially expressed genes was examined throughout expression 
datasets generated from different plants organs (all available 
material, except data generated from cell culture) and sugar 
treatments (AT-0006; 0014; 0015; 0056; 0133; 0199; 0209; 0281; 
0639; 0650.) using the Genevestigator V3 software [18]. With this 
software, hierarchical clustering based on the Pearson correlation 
of the expression profiles of the selected genes was conducted on 
the extracted data. 

Microarray data have been deposited in NCBI's Gene 
Expression Omnibus and are accessible through GEO Series 
accession number GSE46510: (http://www.ncbi.nlm.nih.gov/ 
geo/query/acc.cgi?acc = GSE46510). 

Carbohydrate analyses 

Frozen cell material was ground in liquid N2. The soluble sugars 
sucrose, glucose and fructose were analyzed in ethanol-water 
extracts. Sugars were extracted in 80% ethanol containing 4 mil 
Hepes, pH 7.5, at 80°C for 15 min. After 10 min centrifugation at 
full speed the supernatant was collected, stored on ice and the 
extraction was repeated. Subsequendy the extraction was repeated 
with 50% ethanol containing 4 mM Hepes, pH 7.5. All superna- 
tants were combined and subsequently analyzed for sugar contents 
using a NADP-coupled enzyme assay. The remaining pellets were 
resuspended in 200 mM citrate (pH 4.6), incubated for 20 min at 
95°C, cooled down and digested over night at 55°C by 
amyloglucosidase (Roche). After a subsequent centrifugation at 
full speed for 10 min, the supernatant was transferred and 
analyzed for sugar contents using a NADP-coupled enzyme assay. 
In brief, first all glucose was consumed by the enzymes glucose-6- 
phosphate-dehydrogenase (Roche) and hexokinase (Roche) in the 
reaction buffer (100 mM imidazol, 1.5 mM MgC12; pH 6.9, 
0.5 mM NADP; 1.25 mM ATP) and the OD 344 was monitored 
until a plateau. Subsequendy all fructose was consumed by the 
addition of the enzyme phosphoglucoisomerase (Roche), which 
could be monitored by the increase of the OD 344 . After reaching a 
plateau, finally all sucrose was consumed through the addition of 
the enzyme invertase (Roche). OD 344 was monitored until a 
plateau was reached. Based on the differences between the OD 344 
at the plateau-stages for each sugar, the sugar concentration for 
glucose, fructose and sucrose were calculated in |j,mol hexose/ 
gFW. 

GUS staining 

13-day-old whole seedlings were transferred into GUS-staining 
solution (1 mM X-Glc (Fermentas/Thermo Scientific); 50 mM 



NaHP0 4 , pH 7.0; 1% v/v TritonX-100). The samples were 
incubated overnight at 37 °C. Afterwards the samples were 
washed twice for 2 min in 50% ethanol, twice for 2 min in 90% 
ethanol and finally rinsed and stored in sterile MilliQ;water. The 
samples were analyzed by taking pictures using a Canon EOS 
650D. 

Results and Discussion 

Novel cell system to identify and investigate sugar- 
responsive genes 

In order to identify sugar-responsive genes independently from 

(i) the developmental and cellular heterogeneity of plant tissues, 

(ii) the differential accessibility to sugars and (iii) the presence of 
the studied sugar as energy source, a simplified biological system 
was adapted using habituated stem-cell-like suspension cultures 
from A. thaliana [17]. The system facilitates a homogeneous and 
direct response towards externally applied primary sugars, while 
excluding tissue- or development-specific signals that could affect 
the sugar-dependent gene expression. Contrary to most cell 
cultures, the habituated cells have the capacity to grow and divide 
in a hormone-free medium and possess the ability to differentiate if 
provided with specific hormones [17]. This makes them an ideal 
model to study cellular response to a specific factor. 

To effectively decipher between sugars used as energy source or 
as signals, cell suspensions were screened for growth on Murashige 
and Skoog (MS) media supplemented with 87 mM of different 
sugars. The pentose xylose (Xyl) was identified as a usable carbon 
source, along with Sue, Glc, Fru, trehalose, raffinose, maltose and 
lactose, providing sufficient energy in the culture media for 
sustained growth over a 7 day growth period (Fig. S1A). In 
plants, Xyl can be converted to Glc via the oxidative pentose- 
phosphate pathway [19] and the Xyl-based growth was earlier 
reported both for cultured plant cells [20] and for Arabidopsis 
seedlings kept on Xyl-supplemented growth media [21]. In our 
system, the growth of cells conditioned to grow on Xyl-containing 
media was about 2.5-fold slower than those conditioned to grow 
on Sue (Fig. 1A), but cell length and general appearance of the 
cells under those two conditions were similar (Fig. S1B,C). Along 
with the decreased growth, the Xyl-grown cells had reduced 
contents of internal soluble sugars, with Glc, Fru and Sue levels 
decreased by 20-, 4- and 9-fold, respectively, when compared to 
Suc-grown cells (Fig. 1B-D). Thus, under long-term low-carbon 
conditions, the cells adapted by establishing a decreased growth 
rate, similar to a whole plant response when kept under prolonged 
low carbon regime [3]. 

Deciphering the role of a sugar as a signaling molecule from its 
role as energy-source is one of the main problems hampering 
studies on sugar-mediated regulation of gene expression. Another 
difficulty is a frequently observed rapid conversion of an externally 
supplied sugar, especially Sue, to other sugars [12,22]. We have 
addressed those difficulties by using the Xyl-grown cells, which 
were exposed to low or very low sugar concentrations (1 and 
5 mM) for a short time (up to 3 h). The results suggested that 
treatment with as little as 1 mM Sue for 1 h did not increase the 
specific sugar pools (Fig. S2). Those conditions were also likely to 
limit the use of sugars as energy source. Next, we assessed the 
impact of externally added sugars on expression of BT2 (encoding 
a BTB- and TAZ-domain containing protein), whose expression 
was reported to be regulated by sugar and diurnal rhythms in 
intact plants, with maximum expression in the dark or in carbon- 
starved seedlings [7,23]. Using qPCR analysis, BT2 expression in 
Xyl-adapted cells was sensitive to exogenous Sue application 
(Fig. 2A), with Sue concentration as low as 1 mM being sufficient 
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Figure 1. Effects of Xyl and Sue, as C-sources in the growth media, on packed cell volume (PCV) and intracellular soluble sugar 
concentration in A. thaliana cell culture. (A) Cells-grown in 1x MS medium supplemented with Xyl or Sue. (B-D) Internal concentrations of Glc, Fru 
and Sue after 7 d of growth on Xyl or Sue media. Error bars represent the standard deviation calculated from 3 biological repeats. 
doi:1 0.1 371 /journal.pone.01 0031 2.g001 



to significantly decrease BT2 expression within 1 h (Fig. 2A, box). 
A similar down-regulation was also observed after 1 h treatments 
with either 1 mM Glc or Fru (Fig. 2B), suggesting that BT2 
responds both to Sue and hexoses, and that growth on Xyl does 
not impair BT2 responses to primary sugars. The rapid response of 
BT2 expression to a very low sugar concentration prompted us to 
seek other genes which may rapidly respond to the same sugars 
under the conditions established for the cell culture. 

Functional genomics identification of genes rapidly 
regulated by low sugar concentrations 

To identify genes rapidly responding to low concentrations of 
specific sugars, expression analyses using the Affymetrix ATH1 
microarrays were performed on Xyl-based cell-culture samples 
treated for 1 h with 1 mM of either Glc, Fru, Sue or Xyl in 3 
independent replicates for each treatment. Gene expression was 
calculated as the ratio of the expression value between the 
treatment with a distinct sugar and the values from Xyl-treated 
control samples. Xyl was used as a control, rather than e.g. 
mannitol, since it is penetrable and serves as energy source 
(Fig. 1), similar to Glc, Fru and Sue. Selection of significandy 
expressed genes was performed using stringent criteria (RMA 
normalization and the threshold p-values <0.002) which identified 



a total of 290 sugar responsive genes (Dataset SI) with 
expression ratios ranging between 1.52 and 0.71. The genes could 
be divided into 7 sub-groups, corresponding to those responsive to 
Suc/Fru/Glc, Fru/Glc, Suc/Glc, Suc/Fru, Sue, Glc, and Fru, 
respectively (Fig. 3A). Hierarchical clustering, based on complete 
Pearson correlation of the expression profiles of the identified 
genes, revealed distinct clusters exhibiting both a specific direction 
of expression patterns and responsiveness exclusively to distinct 
sugar molecules (Fig. 3B). Out of the 290 sugar-responsive genes, 
168 (58%) were down-regulated and 122 (42%) - up-regulated 
(Fig. 3A). A similar bias toward down-regulation was frequendy 
observed in studies on auxin- and other hormone-responsive genes 
(e.g. in [24]). From the sensing and signaling perspective, 
suppression of gene expression may represent a faster response 
mechanism than up-regulation of the expression. 

The comparison of the 290 candidate genes with three 
publically available datasets focusing on the transcriptional 
response to (i) sugars and circadian regulation during the diurnal 
cycle [25], (ii) changes in endogenous carbon status [4], and (iii) 
the supplement of 15 mM Sue to C-starved seedlings [7] revealed 
a substantial overlap of 20% only with genes identified upon Suc- 
supplement (Dataset SI). The overlapping genes were members 
of all 7 sub-groups (Fig. S3), suggesting that the in planta response 
to Suc-treatment might have been partly induced by its hydrolytic 




Figure 2. A rapid response of BT2 expression to low doses of sugars in 7-d-old Xyl-grown A. thaliana cell culture. (A) BT2 expression in 
the cell culture supplemented with 0 (dash-dot line), 1 mM (straight) and 5 mM (dashed) Sue for 0, 5, 15, 30 min and 1 and 3 h (n = 3). The box 
represents data for 1 h treatment with Sue. (B) Effects of 1 mM Sue, Glc Fru and Xyl on 672 expression. Mannitol (Man) at 1 mM was a control. 
Significance: f-test; *c< = 0.05, **a = 0.01, n = 7. 
doi:1 0.1 371 /journal.pone.01 0031 2.g002 
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Figure 3. Sue, Glc and Fru induce sugar species-specific changes in global gene expression in A. thaliana cell culture. (A) Treatment of 
7-d-old Xyl-grown cell culture for 1 h with 1 mM Sue, Glc and Fru induced differential expression of 290 genes, which were divided into 7 subgroups, 
representing their responsiveness to a specific set of sugars. (B) Hierarchical clustering of the expression profiles based on Pearson correlation of the 
290 genes resulted in 8 clusters with 4 clusters each containing up- or down-regulated genes. The single clusters represented expression profiles 
responsive to a specific set of sugar species. Significance: f-test; * a = 0.01; n = 13 to 59, depending on the cluster. (C) qPCR expression analysis of 
selected genes in response to a specific sugar (at 1 mM for 1 h). Mannitol (Man) at 1 mM was a control. The numbers in parentheses indicate the 
correlation between the gene expression profiles measured by the microarray and qPCR analysis. Significance: f-test; *a = 0.05, **a = 0.01, n = 8. See 
also Fig. S7 for examples of other genes. 
doi:1 0.1 371 /journal.pone.01 0031 2.g003 



products Glc and Fru. Similar to our microarray data, the majority 
of those overlapping genes (47 out of 60) were down-regulated. 
Further analyses of publically available datasets on intact plants 
treated with Sue or Glc or exposed to extended dark conditions 
again revealed some overlap for each of the dataset with the 
candidate genes identified in the present study (Fig. S4). About 
50% of the 290 genes had a 2-fold change in expression in at least 
one of the conditions analyzed. Expression profiles of the identified 
genes were distributed in a plethora of tissues (Fig. S5), with no 
apparent bias toward any tissue. Thus, the expression of sugar 
responsive genes identified in the present study using an initially 
root-derived heterotrophic cell culture was not restricted to roots, 
heterotrophic background, but encompassed both non-photosyn- 
thetic and photosynthetic tissues. 

Among the 290 candidate genes selected by microarrays, 218 
coded for proteins with known function (Fig. S6) . Several of the 
identified sugar-responsive genes, e.g. those involved in endomem- 
brane trafficking, cell wall synthesis, protein degradation, micro- 
tubule-associated, etc., were related to the cell division process. 
Signaling-related sugar-responsive genes included TOR, coding for 
a key integrator of energy status in eukaryotes [26,27] and b^JP63, 
coding for a transcription factor integrating Glc and ABA signals 
in plants ([28]. TOR, among other functions, regulates Glc- 
dependent repression of B T2 [27], a gene also found as sugar- 
responsive in our study (Fig. 2). 

Sugar responsiveness of selected genes identified by microarrays 
was further verified using qPCR gene expression analysis (Fig. 3C, 



Fig. S7). In those analyses, mannitol was used as a control to 
assess possible effects of Xyl on gene expression. Generally, Xyl 
had only small effect on gene expression, with the exception of 
ERF104 and TPS9; however, the direction of this effect was the 
same as for other sugars tested, suggesting that carbon from Xyl 
can enter a common metabolic pathway for other sugars [19]. For 
effects of Glc, Fru and Sue, the qPCR analyses confirmed that the 
expression of the genes was sensitive to either one, two or three 
sugars, indicating that the microarray data were reliable as to the 
identification of sugar-responsive genes. However, in several cases, 
especially for genes responding to a single sugar on microarrays, 
the qPCR analysis revealed that the genes were also sensitive to 
one or two other sugars. This was most likely due to higher 
sensitivity of qPCR when compared to microarray analyses. 
Indeed, the qPCR-quantified expression of 14 genes correlated 
linearly to the microarray-determined expression with an average 
R of 0.82±0.20 and an overall gene expression ratio ranging 
from 0.81 to 1.34 (for microarrays) and 0.13 to 4.12 (for qPCR). 
Thus, the identification of sugar-responsive genes using stringent 
selection criteria for microarrays was fully confirmed using the 
more sensitive method. Overall, the qPCR analysis allowed for the 
identification of genes responding to either all applied sugars, or 
hexoses only, or Sue and Fru, or Fru only (Fig. 3C, Fig. S7). 

Genes responsive to all three sugars included b/JP63 
{At5g28770), TOR (Atlg50030), At5g22920 (encoding putative 
RING-type Zinc finger protein), TPS9 [Atlg23870, class-II 
trehalose-P synthase), XTH30 {Atlg32170, xyloglucan endotrans- 
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glucosylase-related protein), MGD2 (At5g20410, type-B monoga- 
lactosyl-diacyl-glycerol synthase,), and ERF104 (At5g61600, ethyl- 
ene responsive factor), the latter also Xyl-responsive (Fig. 3C, 
Fig. S7). Some of those genes were earlier reported as sugar 
responsive in studies on whole plants, e.g. expression of At5g22920 
responded to Sue application in Arabidopsis seedlings [7] and, along 
with TPS9, to sugar availability during diurnal cycle in intact 
plants [4,29]. The qPCR-identified genes responding only to Glc 
and Fru, but not Sue, included At2g22080 (encoding unknown 
protein), URH1 (At2g36310, purine hydrolyzing protein) and PRKR 
(At3g49160, pyruvate kinase related) (Fig. 3C, Fig. S7). Both 
URH1 and At2g220S0 have not been previously identified as sugar- 
responsive, whereas PRKR was shown to respond to diurnal 
changes in A. thaliam wt and pgml mutant, and displaying lowest 
expression during high C-availablity [29] . Finally, based on qPCR 
analyses, one gene [At3g25400, encoding unknown protein) was 
responsive to Fru and Sue, but not Glc (Fig. S7), and one gene 
{At3g57540, remorin-like) responded to Fru only (Fig. 3C). 

The qPCR approach allowed to verify the effects of specific 
sugars on expression of selected candidate genes, when compared 
to the microarrays. For all genes regulated by more than one 
sugar, as found both by qPCR and microarrays, the effects of 
various sugars were always in the same direction of expression (up- 
or down-regulation), although the magnitude of response could 
differ (Fig. 3C, Fig. S7). This suggested that different distinct 
sugar signals for those genes may converge at the transcriptional 
level. Alternatively, the results may also imply rapid sugar inter- 
conversions/metabolism, resulting in indirect effects. In addition, 
the results could not differentiate between possible spatial sites of 
sugar perception for each gene, i.e. whether outside or inside the 
cell, nor provide a clue about any signaling pathway transducing 
the sugar signal for a given gene. 

Dissecting sugar signals and perception sites with sugar 
analogues 

Functional analysis of the signalling by distinct sugars was 
performed to verify whether gene expression occurred indepen- 
dently from the impact of the sugar molecule as an energy source, 
and/ or resulted from its rapid inter-conversion into another sugar. 
Non-metabolizable sugar analogues for Glc and Sue were supplied 
at 1 mM for 1 h to Xyl-adapted cell suspensions. The Glc 
analogues included L-Glc, which cannot be efficiently transported 
into the cell, thus allowing for the identification of Glc signaling 
perceived outside of the cell [5]; and 2-deoxyglucose (2Dog) which 
is phosphorylated to a dead-end product by hexokinase (HXK), 
thus allowing for identification of genes which are regulated via a 
HXK-dependent pathway [1,5]. The Sue analogues included 
palatinose (Pal) and turanose (Tur), with only Tur able to enter the 
plant cells [30,31]. These compounds were therefore used to 
dissect gene responsiveness to a specific sugar and in relation to its 
site of perception (outside and/or inside the cells). As no effective 
Fru analogues have been reported, an absence of gene response to 
all the above-mentioned analogues could suggest a Fru-dependent 
gene expression. In studies with sugar analogues, mannitol was 
used as a control since, similarly to the analogues, it is non- 
metabolizable. 

Among the candidate genes demonstrated by qPCR to respond 
to all 3 sugars, b^IP63 and MGD2 were responsive only to 2Dog; 
XTH30 - to L-Glc; TPS9 and TOR - to Tur or Pal, respectively; 
At5g22920 - to 2Dog and Tur; BT2 - to L-Glc and Tur; and 
ERF104 - to none of the analogues tested (Fig. 4, Fig. S8A). The 
results imply that expression of b£IP63, MGD2 and At5g22920 
involves the action of HXK, while expression of XTH30 and BT2 
requires an external Glc (both genes) and external Sue (BT2 only) 



signal. Expression of At5g22920 was additionally regulated through 
a Sue-specific pathway, which requires the uptake of Sue, as 
shown by the Tur effect (Fig. 4). A similar Suc-pathway for 
regulation of gene expression can be proposed for TPS9, but not 
TOR expression, as the latter was affected by extracellular Sue 
(based on Pal effects) (Fig. S8A). The fact that the analogues had 
no effect on expression of ERF104 (Fig. S8A) suggests a putative 
signaling by Fru or that the gene is responsive to a downstream 
metabolite arisen from sugar metabolism. 

For genes found by qPCR to respond only to Glc and Fru, but 
not Sue, sugar analogues had no effect on expression of URH1 
(Fig. 4), suggesting that the actual signaling molecule is Fru or 
another signal, which could not be identified with our experi- 
mental setup. A change in expression of At2g22080 and PRKR was 
induced by 2Dog (Fig. S8A), consistent with regulation through a 
HXK-dependent process. Unexpectedly, Pal affected PRKR 
expression (Fig. S8A), in contrast to the lack of Sue effect on 
this gene (Fig. S7) . This may represent an example of a distinctive 
signal induction derived from the non-metabolizable sugar 
analogue (Pal), as reported earlier [32]. Expression of At3g25400, 
which was responsive to Fru and Sue, but not Glc (Fig. S7), was 
affected by L-Glc and 2Dog (Fig. S8A), suggesting that the Sue 
effect was indirect, resulting from its hydrolysis to Glc and Fru, but 
the true nature of the signal is unclear. Finally, At3g57540 
expression, which responded to Fru only (Fig. 3C), was not 
affected by the analogues applied (Fig. 4), suggesting a Fru- 
dependent pathway, similarly to URH1. For all genes where Fru, 
Glc and 2Dog affected the expression, Fru effects were likely 
signaled via a Fru-specific pathway, which is distinct from HXK- 
mediated mechanism, but feeds into the same down-stream 
regulation pathway [8,33]. 

In addition to studies with the above-mentioned sugar 
analogues, we have also tested effects of 3-O-methylglucose 
(30MG), a D-Glc analogue, which is transported into plant cells 
but is metabolized very slowly [34]. This compound, although at 
first frequently used to identify genes responding to Glc via a 
HXK-independent pathway (e.g. in [5]), was later proposed to be 
not perceived as a sugar signal by plants [34] and was ineffective in 
modulating the expression of over 200 Glc-responsive genes in 
Arabidopsis [35]. In our experimental setup, 30MG at 1 mM had 
significant effects on expression of MGD2, TOR and PRKR (Fig. 
S8B); however, the direction of expression induced by 30MG was 
in all three cases opposite to the direction of sugar-induced 
expression for a given gene. Thus, strikingly, 30MG behaved as a 
classical antagonist of Glc effect, suggesting that it interferes with a 
Glc perception/transduction pathway. It is unclear, however, 
whether those results are relevant to sugar signaling, given the 
controversy surrounding the effectiveness of 30MG in truly 
reflecting a sugar signal in plants. 

Studies with sugar analogues have demonstrated that the cells 
had the ability to effectively inter-convert exogenously added sugar 
molecules, and this concerned not only the inter-conversion of Sue 
to/from Glc and Fru, but also between the two hexoses. The inter- 
conversions occurred even though the sugars were applied at low 
concentration (1 mM) and for short time (1 h), i.e. conditions 
which apparendy had little or no significant effect on internal 
sugar pools (Fig. S2). This suggested that the active quantity of 
the exogenously applied sugar required to induce a change in gene 
expression was lower than 1 mM. Nevertheless, using a combi- 
nation of treatments with a given sugar and a given sugar 
analogue, we were able to unequivocally pinpoint the true nature 
of the signalling sugar for a given gene tested. The genes 
responded rapidly to low concentrations of a given sugar and can, 
thus, be considered as part of an early response mechanism to the 
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Figure 4. Effects of sugar analogues on expression of selected sugar-responsive genes in 7-d-old Xyl-grown A.thaliana cell culture. 

The analogues were added at 1 mM for 1 h. Man (mannitol), L-glc (L-glucose), 2Dog (2-deoxyglucose), Pal (palatinose), Tur (turanose). Significance: f- 
test; *a = 0.05, **cx = 0.01, n = 5. See also Fig. S8A for examples of other genes. 
doi:1 0.1 371 /journal.pone.01 0031 2.g004 



sugar availability for cells otherwise adapted to grow under long- 
term low carbon conditions. 

Altogether, the functional dissection with sugar analogues 
allowed to identify genes that are specifically regulated not only 
by one sugar species (e.g. Fru, At3g57540), but also by two (Glc and 
Sue for At5g22920) or, for some genes, perhaps even by each of the 
three sugars, given that we could not unequivocally prove Fru- 
specific effects. The mechanisms of a given sugar signaling may 
also differ, with Sue or Glc sensed outside or inside of the cell, with 
some of the Glc effects signaled via the HXK-dependent or HXK- 
independent pathways. In Table S2, we have summarized the 
results of sugar- and sugar analogue-dependent expression of 
selected genes and have listed the identified sugar-specific gene 
responses and sites of sugar sensing for each gene. 

Expression of sugar-responsive genes in planta 

To assess the biological relevancy of the use of habituated cell 
culture in identifying sugar-responsive genes, expression of selected 
genes was studied first in transgenic plants transformed with 
promoter-GUS reporter. GUS-expression analyses for b%IP63, 
AtRS2.ll, URH1 and TOR revealed the cell/tissue/organ-specific 
expression patterns of each gene: b^IP63 - in roots, cotyledon 
meristems and leaf vasculature, TOR - in meristems of both apical 
shoot and side roots, MRS2.11 - mainly in leaf mesophyll cells, and 
URH1 - in roots only (Fig. 5). A similar distribution of GUS- 
expression for those genes was reported earlier [36-39]. The 
expression of the genes in whole seedlings was not restricted to one 
specific tissue- or cell-type, but represented a variety of tissues/ 
cells, both autotrophic and heterotrophic. This was supported by a 
more detailed Genevestigator-assisted analysis of tissue-dependent 
expression which included - besides b^JP63, MRS2.11, URH1 and 
TOR - also other genes that were selected for qPCR studies (Fig. 
S5B). Thus, the identification of sugar-responsive genes in cell 
culture was not biased toward genes preferentially expressed in 
any specific type of cells/tissues. 

In subsequent studies, using qPCR, we assessed whether the 
genes identified with the cell culture system were also sugar- 
responsive in intact Arabidopsis plants which experienced changes 
in internal sugar concentrations. When wt plants and the starch- 
deficient pgml mutant, both grown in the 16 h light/8 h darkness 
photoperiod, were exposed to a 6 h extended night conditions, the 
endogenous soluble sugar contents decreased in comparison to the 
end of light conditions (Fig. 6A). Changes in sugar contents were 
accompanied by changes in expression of sugar-responsive genes; 
this was observed in both shoots and roots for both wt and pgml 
seedlings (Fig. 6B, Fig. S9). Since the pgml has little or no starch 




Figure 5. The expression of selected sugar-responsive genes in 
planta is not restricted to a specific tissue or cell type, as found 
by GUS expression analyses of four promoter-reporter lines 
for bZIP63 (A), TOR (B), MRS2. 1 1 (C) and URH1 (D) in 1 3-d-old 
Arabidopsis seedlings. 
doi:10.1 371/journal.pone.01 0031 2.g005 

(Fig. 6A), the similar expression patterns of sugar-responsive 
genes in wt and pgml plants suggested that their expression was 
independent of the formation of starch or the availability of starch 
as backup energy source. Generally, we observed little or no 
difference in the contents of soluble sugars between the pgml 
mutant and wt plants (Fig. 6A). Even though these results are 
contrary to observations published earlier, it has been highlighted 
that the ratio of soluble sugar content in the pgml versus wt plants 
is dependent on growth conditions, especially the length of the 
photoperiod (reviewed in [40]). In contrast to other studies, the 
seedlings were germinated and grown in vitro under long-day 
conditions (16 h light/8 h dark regime) which resulted in a 
comparable phenotype of pgml and wt. In addition, the samples 
were harvested at an early developmental stage (after 13 days), a 
time point where the difference in the soluble sugar contents 
between pgml and wt is less established. 

For a given gene, the direction of the change of expression 
(increase or decrease) was as observed in the sugar-treated cell 
culture. Thus, BT2, bZIP63, At5g22920, and TPS9 were signif- 
icantly up-regulated by dark conditions (low sugar content), 
whereas MGD2 and MRS2. 1 1 were up-regulated in the light (high 
sugar content) (Fig. 6B, Fig. S9). Other tested genes, with the 
possible exception of URH1, followed the same trend, even though 
the results were not always statistically significant. Similar to the 
GUS-expression patterns observed in transgenic plants (Fig. 5), 
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Figure 6. Endogenous sugar contents (A) and the expression of cell-culture-selected sugar-responsive genes (B) respond to light/ 
dark conditions in intact plants. Shoots and roots from 13-d-old A. thaliana seedlings of wt (light bar) and pgml (dark bar), grown at 16/8 h 
photoperiod, were harvested at the stages of high (end of 16 h day, white box) and low (end of 6 h extended night, black box) carbon availability. 
Significance: f-test; *a = 0.05, **a = 0.01, n = 3. See also Fig. S9 for examples of other genes. 
doi:1 0.1 371 /journal.pone.01 0031 2.g006 



the highest MRS1 expression was in shoots, while URH1 was 
mainly expressed in roots. Both b£IP63 and TOR were expressed 
about equally in shoots and roots (Fig. 6B, Fig. S9). 

Overall, the results strongly suggest that the photoperiod- 
mediated changes in endogenous sugar contents in planta play a 
significant, if not decisive, role in regulating expression of the 
selected genes. The identities of effective sugars could not be 
inferred from those data, as contents of all sugars changed in a 
similar way during the light/dark conditions (Fig. 6A), and we 
could not control sugar contents within each specific cell-type 
within the organ studied (shoot vs. root). However, for the first 
time, we can rationally infer the true identities of the effective 
sugars and specific sugar-signaling, based on the data obtained 
from the habituated cell cultures. 

Habituated cells as a model system to study sugar 
signalling 

In plants, each cell type either produces its own sugars and/ or 
imports them from other cells. As sugars are the main energy 
source and provide the backbone for all carbon-containing 
compounds, the proper functioning of the whole plant depends 
on a highly coordinated network of sugar molecules being 
produced and metabolized, exported and imported, both between 
adjacent cells and via long-distance transport through the phloem. 
This adds to the complexity of studies on sugar-signaling in whole 
plants, where each cell type or tissue may face different sugar 
concentrations and may have distinct cell/tissue-specific mecha- 
nisms modulating expression of sugar-responsive genes. Another 
difficulty concerns deciphering between the direct and indirect 
effects of sugars, as an increase in sugar content stimulates both 
auxin and ABA synthesis [28,41], whereas sugar starvation inhibits 
gibberellin production [42]. In intact plants, studies on the impact 
of sugars and hormones in modulating gene expression have to 



take into account the respective weight/impact of their cross-talk 
in each cell/ tissue-type within the organism. These limitations do 
not apply when using the habituated cell culture system, which is 
moreover highly versatile with respect to requirements for carbon 
source (Fig. SI). 

Sugar-responsiveness of a given gene depends on availability of 
specific sugar signal(s), but it may also depend on developmental 
status of different cells/tissues. In the present study, sugar- 
responsiveness of the identified genes, although determined within 
the context of a heterotrophic cell culture, was not restricted to a 
distinct cell, tissue or organ in whole plants (Fig. 5, Fig. S5). 
Moreover, expression of those genes was correlated with internal 
changes in sugar concentration in whole plants during diurnal 
cycle, but the magnitude of expression was dependent on a given 
organ (Fig. 6). This strongly suggests that the results obtained 
with the habituated cell culture are of high relevance to studies on 
sugar signaling in whole plants. Aside from other advantages, the 
culture represents an ideal system in which exogenous supply of 
specific sugar(s) and/or morphogen(s) and/or cell differentiation 
status can be controlled by the experimenter. The system can also 
be used to study effects and mechanisms of action of any other 
signaling compounds (e.g. a plant hormone), with or without the 
sugar cross-talk. 

Supporting Information 

Figure SI Effects of various C-sources on growth of A. 
thaliana cell culture (A) and specific effects of Xyl (B) 
and (C) on cell length and general phenotype of the 
cultured cells. 

(DOCX) 
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Figure S2 Effects of Sue concentration and induction 
time on the intracellular contents of soluble sugars in A. 
thaliana cell culture. 

(DOCX) 

Figure S3 The global gene expression response to 
exogenously applied sugars overlapp between A. thali- 
ana cell culture and seedlings. 

(DOCX) 

Figure S4 The effect of the treatment of plant material 
with varying Sue and Glc concentrations on the expres- 
sion of 290 identified genes (A) and of the 14 selected 
genes (B). 

(DOCX) 

Figure S5 Distribution of the expression of 290 identi- 
fied genes (A) and of the 14 selected genes (B) in various 
tissues/organs. 

(DOCX) 

Figure S6 Functional classification of 290 identified 
genes. 

(DOCX) 

Figure S7 Defined sets of sugars trigger significant 
changes in gene expression of selected sugar-responsive 
genes in 7-d-old A. thaliana cell culture. 

(DOCX) 

Figure S8 Effects of sugar analogues on expression of 
selected sugar-responsive genes in 7-d-old Xyl-grown 
A.thaliana cell culture. 

(DOCX) 
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